Theory of momentum resolved tunneling into a short quantum wire 
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Motivated by recent tunneling experiments in the parallel wire geometry, we calculate results 
for momentum resolved tunneling into a short one-dimensional wire, containing a small number of 
electrons. We derive some general theorems about the momentum dependence, and we carry out 
exact calculations for up to N = 4 electrons in the final state, for a system with screened Coulomb 
interactions that models the situation of the experiments. We also investigate the limit of large 
iV using a Luttinger-liquid type analysis. We consider the low-density regime, where the system is 
close to the Wigner crystal limit, and where the energy scale for spin excitations can be much lower 
than for charge excitations, and we consider temperatures intermediate between the relevant spin 
energies and charge excitations, as well as temperatures below both energy scales. 
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I. INTRODUCTION 

Mesoscopic and nanoscale systems have proved to be 
fertile grounds for studying quantum systems where in- 
teractions play a vital role in the physics. In particular, 
superconductivity*^ the Kondo effectj 3 ^^ and Luttinger 
Liquids (SiSuLiii have all been studied experimentally and 
theoretically to varying degrees in these systems of re- 
duced size and dimensionality. Yet, due to the myriad 
of energy scales that confinement and interactions intro- 
duce, there is an ever unfolding landscape of rich physics 
to be explored and understood. Many of these interest- 
ing regimes in mesoscopic and nanoscale systems can only 
be explored experimentally with advances in technology. 
One such advance was the ability to fabricate quantum 
wires of extremely high quality in so-called cleaved edge 
overgrowth samples: 11 Auslaender et aim took this a step 
further to fabricate two parallel quantum wires of very 
high quality which permitted momentum to be approx- 
imately conserved in the tunneling between them. See 

Fig.n 

Since then a series of experiments has been car- 
ried out at the Weizmann Institute by Yacoby and 
collaboratorsS&iS on such cleaved edge overgrowth sam- 
ples fabricated on a GaAs substrate. The momentum 
transfer in the tunneling was controlled by application 
of a magnetic field perpendicular to the plane containing 
the two wires; the energy-dependence was explored by 
varying the voltage between the wires; and the electron 
density could be varied by application of voltage to a top 
gate. With proper interpretation, these tunneling mea- 
surements give very detailed information about the quan- 
tum wires. For example, one can measure directly the 
Fermi wavevectors associated with the occupied modes 
in the two wires; one sees diffraction fringes due to the fi- 
nite length of the shorter of the two wires, from which one 
can deduce the shape of the gate potential that confines 
the electrons at the ends of the shorter wireji 3 - and one 
can see evidence of the strong effects of electron-electron 
interactions in one-dimensional systems, manifest in dis- 



tinct propagation velocities for spin and charge, as pre- 
dicted by Luttinger-liquid theoryiVWWV 8 Earlier 
experiments on the conductance of finite quantum wires 
motivated a number of theories in which the finite size 
effects were explored in the context of Luttinger Liquid 
theory, 19 ! 20 . 21 ! 22 . 23 

In a recent extension of the tunneling experiments, 
Steinberg et alM^ have constructed a pair of wires in 
which a central portion, of order 2 /zm in length, is cov- 
ered by a gate and can be depleted of electrons relative to 
the rest of the wire. Measurements revealed a dramatic 
transition, when the density of electrons in the center of 
the upper wire was reduced below a critical value. The 
data give strong evidence that in this regime, the elec- 
trons remaining in the depleted wire segment are sepa- 
rated from the rest of the wire by barriers at the two ends, 
so that the number N of electrons in the wire segment 
becomes quantized in integer values. For example, the 
tunneling conductance showed a series of sharp maxima, 
as the gate voltage was varied, analogous to Coulomb 
blockade peaks, which were associated with transitions 
between N and N — 1 electrons in the ground state of 
the wire segment. The positions of the Coulomb block- 
ade peaks are independent of the applied magnetic field, 
but there is variation of the peak height with field, which 
enables one to measure the momentum distribution of 
the many-body quantum state at the Fermi level. The 
momentum distribution is quite different from what is ob- 
served at electron densities above the critical density, and 
it was suggested that in the low-density regime, the elec- 
trons are rather strongly localized by electron-electron 
interactions as well as by the barriers at the end of the 
wire segment. It should be mentioned that this regime of 
localized states was observed not only near depletion of 
the lowest mode of the upper wire, but also as the elec- 
tron density in the second and third mode was reduced 
below a critical value. 

At the present time, we do not have a good under- 
standing of why a potential barrier should arise sponta- 
neously between the low density central segment and the 
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FIG. 1: Schematic geometry of electron tunneling between 
two parallel quantum wires. Electrons are assumed to tun- 
nel between an infinitely long lower wire and a short upper 
wire of length L. The wires are separated by a center-to- 
center distance d. The upper (lower) wire has a width W u 
(Wi) and an average electron density n u (n;). A magnetic 
field B is applied perpendicular to the plane of the wires 
to allow a field-dependent momentum boost hqs = eBd in 
the tunneling. The energy of the tunneling electrons can be 
changed by adjusting the voltage applied between the two 
wires. In the experiments of Ref. Q, d = 30 nm, W u ~20 nm, 
Wi=30 nm and L — 2 — 10 (j,m. Electron densities can range 
n = — lOOjLtm - . For much of this paper we will assume 
(nzas) -1 <C 1 and (n^as) -1 ^> 1, where as is the Bohr ra- 
dius, so that the upper wire is close to the Wigner crystal 
regime and the lower wire can be treated as non-interacting. 



higher density ends of the upper wire. Barriers of this 
type have been found in Hartree-Fock calculations of a 
one-dimensional model 26 that is a crude representation 
of the geometry of the experiments of Ref. 0E3i as 
well as in earlier LDA calculations of a quantum point 
contact. 27 However, it remains to be seen whether these 
results are robust or are artifacts of the particular models 
and approximations. 

In any case, inspired by the results of Steinberg 
et ai-r^Si we have been led to consider the theory 
of momentum-resolved tunneling into a short quantum 
wire, containing a small number of electrons. We con- 
sider tunneling from an infinite lower wire into a short 
upper wire, under conditions where the energy levels of 
the upper wire are discrete, and the temperature is small, 
at least compared to the Coulomb blockade energy and 
the Fermi energy in the upper wire. We will be particu- 
larly interested in the situation where the density in the 
upper wire is small, so that when two electrons come close 
together, the electron-electron repulsion is strong com- 
pared to the average kinetic energy, and it is difficult for 
two electrons to change places. At low densities, it can 
be helpful to think about the one-dimensional electron 
system as a kind of fluctuating Wigner crystal^^S^ with 
a Heisenberg antiferromagnetic exchange coupling32' 3 i2 2 i 
J between successive spins that is small compared to the 
energy for short-wavelength charge fluctuations or the 
original Fermi energy. Long-range interactions are ex- 
pected to enhance the Wigner crystal-like correlations, 33 
and further reduce J. 

For a finite low-density system, we can go one step 
further and consider a situation where J is so small that 
the temperature T can be larger than J/ks but still small 
compared to the lowest energy charge excitation, whose 
energy is proportional to the charge-propagation veloc- 



ity divided by the length of the wire. This situation, 
which we denote the "free-spin regime," will be consid- 
ered in this paper, along with the strict low-temperature 
regime (the "ground-state tunneling regime" ) where fcgT 
is small compared to all excitation energies. We shall 
also discuss a situation kgT 3> NJ, in which fc^T is 
larger than all spin excitation energies, including multi- 
ple spin excitations, which we shall refer to as the "ex- 
treme free-spin regime." We shall further distinguish be- 
tween a "true Wigner crystal" and a "fluctuating Wigner 
crystal" depending on whether the root-mean-square dis- 
placement of an electron in the center of a finite wire, 
due to quantum fluctuations, is smaller or larger than 
the mean electron spacing. 

In our discussion, we shall treat the infinite lower wire 
as a collection of non-interacting electrons. See Fig. ^ 
Interactions between the two wires are considered only 
in so far as the lower wire may act as an ideal conductor 
that contributes to a frequency-independent screening of 
the Coulomb interaction between electrons in the upper 
wire. The high-density end-sections of the upper wire are 
taken into account only as weakly coupled leads, which 
set the chemical potential in the confined central part of 
the wire. 

We shall concentrate on tunneling in the limit of zero 
bias voltage, and shall examine both the momentum de- 
pendence (i.e., the dependence on applied magnetic field) 
and the overall magnitude of the tunneling conductivity, 
when the gate voltage is adjusted to give degeneracy be- 
tween the ground states of the upper wire for N and 
N — 1 particles. Much of our discussion will focus on 
cases where TV is quite small, but we shall also consider 
the behavior as N becomes large. 

We assume that the coupling of the upper wire to its 
leads is weak enough that the resulting level broadening 
is small compared to the thermal energy ksT. However, 
we assume that the conductance between the upper and 
lower wires is weaker still, so that the resistance is dom- 
inated by the tunneling conductance, and so we can cal- 
culate the current using Fermi's golden rule, in which the 
tunneling matrix element enters only to second order, as 
a prefactor to the overall conductivity. 

This paper is organized in the following way. In Sec.lTTl 
we introduce the basic model and notation that we use 
throughout the paper. In Sec. IIIII we discuss the impor- 
tant special case of non-interacting electrons in the upper 
wire. In Sec. II VI we prove some general theorems about 
the momentum structure of the tunneling for small N. 
In Sec^we present a numerical study of the double wire 
geometry, including interwire screening effects, for small 
N based on exact diagonalization. In Sec. I VII we discuss 
the general features expected for tunneling in the limit 
of large N, focusing on the free-spin regime, where the 
interactions are very strong and the temperature is large 
compared to the spin exchange energies. In Sec. IVIll we 
compare our theoretical findings with the measurements 
by Steinberg et al.^^ before summarizing the paper in 

secEnn 
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II. MODEL AND NOTATION 

The Hamiltonian for electrons in the upper wire con- 
sists of a kinetic energy, a one-body confining potential 
U(x), and a two body interaction v(x — x'), where x is 
the distance along the wire. We shall discuss some theo- 
rems, below, which arc independent of the detailed form 
of U and v. In our more detailed calculations, however, 
we shall concentrate on the case where U contains infi- 
nite vertical walls at x = and x — L, together with a 
potential arising from a uniform positive background in 
the interval < x < L. The interaction v will have the 
form of a Coulomb potential at intermediate distances, 
screened at long distances by the parallel lower wire, and 
cut off at short distances because of the finite thickness 
of the upper wire. With these choices, one can consider 
a limit N — > oo,L — * oo with N/L fixed, and the elec- 
trons will be uniformly spread along the wire. We shall 
primarily be concerned, however, with systems where N 
is not very large. 

Let ty^ and Vl/^ -1 represent many body states of the 
upper wire with N and N— 1 electrons, respectively, with 
energies and E^~ l . We consider here the case where 
the voltage V between the upper and lower wire is small 
compared to fcgT, and the tunneling conductance may 
be written as G = G+ + G_ in terms of 



G± = CB{k±) , 



where 



TO = £l(*«|cLl* 



N\J l,T/W-l\|2 



Wo 



(1) 

(2) 



In the limit where fc^T is small compared to any ex- 
citation energies of the system, only the ground states 
contribute to the sum in J2J, and one obtains a large 
value of G only if the gate potential is adjusted so that 
E — E 1 ^^ 1 — fi is close to zero, i.e., of order ksT or 
smaller. If E N — E N ~ X = fj, precisely, we have 



G, 



Afi 



X 2 (3vL 



1 



9n + 9n-i 



B(k d 



where 



B(k) 



E 

ot-ycr 



JV-1\ |2 



(7) 



(8) 



and <7jv and gN-i are the degeneracies of the ground 
states V N and f^ -1 . 

For the system under consideration, with spin- 
independent interactions and no spin-orbit effects, eigen- 
states of the Hamiltonian will have a definite quantum 
number for the total spin S. In the absence of a Zeeman 
field, then, the iV-electron ground-state will have a de- 
generacy <7jv = (2Sjy + 1), where Sn is the corresponding 
spin quantum number. The indices a and 7 may then 
be chosen to label states with different values of the spin 
component S z . In the case where there is an applied 
magnetic field with a Zeeman energy large compared to 
ksT, however, the degeneracy will be removed, and we 
need only consider the states with spins aligned parallel 
to the field, S z = S. 

When the Zeeman field is absent, we may use the 
Wigner-Eckart theorem to carry out the sum over the 
spin indices, and write 



/(£a 7 ) 



c a7 



= E, 



N 



E. 



N-l 



-PUN 



N-l 



(3) 



(4) 



(5) 



and cj.^ = / Q L dxipl(x)e lkx / \/L creates a particle in the 
upper wire with spatial wave function e lkx and spin 
(7 = ±1/2, while v is the (total) density of states per 
unit length at the Fermi energy of the lower wire, A 
is the coefficient of the tunneling Hamiltonian (properly 
normalized), \x is the chemical potential, f3 — l/(fcsT), 
/(e) = 1/{1 + exp[/3(e — fi)}} is the Fermi function, and 
Zn is the partition function of the upper wire when it 
contains precisely N electrons. (We neglect the proba- 
bility of occupancies other than N or N — 1.) k± are 
given by the lower wire wavevectors at the Fermi energy, 
shifted by the perpendicular magnetic field B: 



±k[ 



eBd/h 



(6) 



B{k) = D\M(k)\ : 



where M is the matrix element 



M(k) = 



N-l\ 



(9) 



(10) 



with the indices for S z chosen as a = Sn and 7 
and a — a — 7. The constant D is given by 

D = 1 + 2max(SW,SAr_i) . 



S 



N-l, 



(11) 



In order to have a non-zero matrix element, we must have 
|iSjv — Sn-i\ = 1/2. This condition is always satisfied for 
the ground states, because for a system obeying the one- 
dimensional Schrodingcr Equation with spin-independent 
interactions, the lowest energy state necessarily has spin 
quantum number S = 0, if N is even, and S — 1/2, if 
N is odd, according to the Lieb-Mattis theoreniiM This 
theorem also tells us that D = 2 for tunneling between 
ground states^ 

It is instructive to define a "quasi-wavefunction" 



(12) 



so that M(k) = J dxe lkx ^ s *(x)/VZ. For non- 
interacting electrons, ^f^ s (x) is the wave function for the 
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last added electron, and its norm is 1. For interacting 
electrons, the norm will be less than 1, due to orthogo- 
nality catastrophe-type effects. 

We remark that the value of G given by equation Q is 
not actually the maximum value of the conductance that 
can be achieved at the given temperature T. A somewhat 
larger conductance can be achieved by shifting the gate 
voltage slightly away from the value where the states for 
N and N — 1 particles have precisely the same energies, 
in order to take advantage of the differing degeneracies 
of the states. One finds that the maximum conductance 
is related to G in Eq. J7J) by 



In either case, one finds that 



Gmax/G = 2( J ? 2 + l)/(l+7 ? ) 2 , 



(13) 



where rj — (g^ / gjv-i) 1 ^ 2 • The maximum occurs when 
E N — E N ~ X — /i = fcsTlnry, so that there is equal proba- 
bility for the upper wire to have N or N— 1 electrons. For 
ground-state tunneling between states with spins and 
1/2, the correction factor for the maximum conductance 
is G max /G = 1.0294. 

If the spin excitation energies are sufficiently low com- 
pared to all other excitation energies, we may consider 
tunneling in the "extreme free-spin" regime, where ksT 
is large compared to the spin energies, and we may ne- 
glect the energy differences between different spin states. 
The degeneracy factors in the expression for C, Eq. JSJ, 
are then given by gpj — 2 N , and the labels a and 7 in 
the equation for B(k), Eq. (J2J), must be summed over all 
of the states in the manifold. 

In the formulas above, we have assumed that the local- 
ized section of the upper wire is separated from its high- 
density leads by barriers sufficiently high that the level 
width r due to coupling to the leads is small compared 
to ksT. In the other limit, where fcgT -C T, the pref- 
actor C for ground-state tunneling will be proportional 
to r _1 rather than T -1 . If eV is larger than both ksT 
and r, then we cannot use the linear conductance for- 
mulas, but the dependence of the current on momentum 
transfer should still be proportional to B(k) at moderate 
voltages. 



III. NON-INTERACTING ELECTRONS 

It is useful to begin by recalling some results for 
non-interacting electrons. We consider particularly the 
case of square-well confinement, where U(x) = for 
< x < L, and 00 otherwise. The one-electron states 
are then sine waves, with spatial dependence 4>n(%) — 
(2/L) 1 / 2 sva{kjjx], where kfj = Ntt/L , and, in the ab- 
sence of a Zeeman field, 



N 



N+l 



(14) 



Int 



with [XJint the integer part of X. For fully polarized 
electrons (or spinless fermions), we have simply N = N. 



M(k) = e lkL/2 



L 



i[(k + k^)L/2] 



N 



x sm[(k- k„)L/2] 
k — kjy 



(15) 



We note that the matrix element vanishes at k — if 
and only if N is even. 




FIG. 2: (Color online.) Momentum resolution of one-particle 
box states. Shown is the modulus squared of Eq. 115H for var- 
ious values of N . Many of the features in the non-interacting 
case remain when interactions are included and the Fourier 
transform of a many-body wavefunction is taken. 

For large N, we see in Fig. |2 that |M(fc)| 2 has sharp 
maxima at k — i/c^y, which become (^-functions in the 
limit N — > 00, L — * 00. Of course, k^ corresponds to the 
Fermi momentum of the infinite system. For large but fi- 
nite N, the (^-functions broaden into diffraction patterns, 
with spacing 8k = 2w/L between successive zeroes, and 
a peak-width of the same order. 

Tserkovnyak et al&H have discussed the form of M(k) 
for large but finite N when the confinement at the end 
of the wires is soft on the scale of the Fermi wavelength. 
In this case they find an asymmetric diffraction pattern, 
which falls off rapidly for |fc| > fcjy, but relatively slowly 
on the side |A:| < fejy. The spacing between zeroes is not 
constant and is slightly larger than 2-njL. 



IV. SOME THEOREMS FOR SMALL SYSTEMS 

Some rigorous theorems are helpful for understanding 
the properties of B(k), particularly at small values of N. 

As a consequence of time-reversal invariance, all wave 
functions may be chosen to be real. From this it follows 
that B(k) — B(—k). Also, the matrix element M(k), 
which determines B{k) in the case of ground state tun- 
neling, satisfies M(k) = M(-k)*. 
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If the confining potential is symmetric about the center 
of the well, the states ^> N and vl/^ -1 may be classified 
by their parity under inversion through the center of the 
well. If the states have the same parity, then M{k) = 
M(—k), and hence M must be real. If the states have 
opposite parity, then M(k) = —M{—k), and M must be 
purely imaginary. In either case, if M(k) changes sign, 
as a function of k, there will be one or more values of 
k where M(k) = 0. Thus, in the ground-state tunneling 
regime, there can be points k where B(k) = 0. We have 
already seen an example in the case of non-interacting 
electrons. 

For the case N = 1, the ground state wave function 
can be written in first quantized notation (in terms of 
position x and spin index a), as 



(16) 



where ^(x) > for all x. (This is true regardless of 
whether or not the confining potential has inversion sym- 
metry.) Thus, for the case N = 1, the matrix element M 
can be written as 



M(k) 



dx—^^/*(x) 



(17) 



Clearly, |M(fc)| has its maximum when k = 0, for N = 1. 

We obtain a similar behavior for ground-state tunnel- 
ing in the case N = 2. The ground state of the two- 
fermion problem, in any number of dimensions, is always 
a spin-singlet, with a spatial wave function that has no 
nodes. Thus, for N = 2 we may write, in first-quantized 
notation, 



(18) 



where Xo = 2~ 1 / 2 (| fj,) — | 11)), while / is real and posi- 
tive, and symmetric under interchange of its arguments. 
We may also write 



4J*T> = 



[e ik ***( Xl ) \H)-e tkx ^(x 2 ) ||T>], 



where ty(x) > is the wave function for the one-particle 
system. The matrix element M(k) is then given by 

M(k) = J dx 1 dx 2 ^-y(x 1 )f(x 1 ,x 2 ) , (20) 

which clearly has its maximum at k = 0. 

The situation is different if we consider tunneling be- 
tween the N = 1 system and the lowest triplet state of 
the N = 2 system. If we choose the spins to be maximally 
aligned in the z-direction, we have 



\*)=g(x l ,x 2 )\ 11), 



(21) 



and 



4 t I*t) 



where g is antisymmetric under interchange of its argu- 
ments, and 



[e lkx ^{x l )-e lkx ^{x 2 )] | TT> , (22) 



M(k) = V2 I dx 1 dx 2 —^^(x 1 )g(x 1 ,x 2 ) ■ (23) 



If the confining potential is symmetric, then g will be 
antisymmetric under inversion, while "J is symmetric. It 
follows that M = at k — 0, and |M(/c)| will have its 
maximum at a non-zero value of k. 

The triplet state will be the actual ground state of the 
N = 2 system if there is an applied Zccman field large 
enough to overcome the exchange splitting between the 
singlet and triplet states. It is also the correct ground 
state for a model of spinless fermions. In the absence of 
a Zeeman field, in the free-spin regime, where the temper- 
ature is larger than the exchange splitting, but smaller 
than all other excitation scales, the quantity B(k) will 
have two contributions from the singlet state and four 
contributions from the triplet. In this case, the function 
B(k) will not be zero at k = 0, but we may have a situ- 
ation with B(k) having a local minimum at k = 0, with 
a pair of maxima at a finite value of |fc|, as in Fig. 1101 
below. 

For N = 3, the many-body ground state for non- 
interacting electrons in symmetric confinement has odd 
parity under inversion, as it has two electrons in the low- 
est even-parity level, and one electron in the second odd- 
parity level. The N = 4 state has even parity. These par- 
ities cannot be altered by a weak electron-electron inter- 
action, and it is likely that they will hold even for rather 
strong interactions. As we have seen, the ground state 
for TV = 2 has even parity. Consequently, for electrons 
in a symmetric confining potential, the ground state tun- 
neling form factor B(k) should vanish at k = for N = 3 
and TV = 4 . 

The feature of a zero amplitude at k = does not 
persist, if the confining potential is not symmetric under 
inversion. For a sufficiently asymmetric potential there 
need not be even a local minimum near k = 0. For ex- 
ample in the case of an asymmetric double- well potential, 
when there is little overlap between the wave functions 
in the two wells, one would expect that B(k) for N = 3 
and N — 4 should be qualitatively similar to the form for 
N = 1 and N = 2. 

The state N = 5 has even parity, at least for weakly 
interacting electrons. Thus the matrix element for tun- 
neling between the ground states with 4 and 5 electrons 
will be non-zero at k = 0. For non- interacting electrons, 
B(k) has a local maximum at k = 0, but has its abso- 
lute maxima at a non-zero value of |fe|, with zeroes of the 
amplitude in between. See Fig. [3 We expect these quali- 
tative features for non-interacting electrons to persist for 
strong interactions in small- TV wires. 

We now turn to exact numerics for small N. 
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V. CALCULATIONS FOR N = 1 TO 4 

We can illustrate some of the ideas and results dis- 
cussed in Sec. IIVI with exact numerical calculations for 
small N which take into account finite size and screening 
effects from the lower wire. 



A. Form of electron interaction 

We assume the interaction potential is a Coulomb in- 
teraction with both a short and a long range cutoff. The 
short range cutoff is achieved in a simple way by modify- 



ing a - potential to 



where z — x — x' is the sep- 



aration of two electrons along the wire and W is a short 
range cutoff (W u for the upper wire), which is roughly the 
width of the quantum wire. The long range cutoff, with 
experiments by Steinberg et alM^ in mind, is achieved 
by putting a second, lower wire with width Wi, assumed 
to be perfectly conducting, parallel to the system under 
study, at center to center distance d. See Fig. ^ This 
approximation is valid, strictly speaking, if the electron 
density in the lower wire is sufficiently high that many 
modes are occupied, while there is only one mode occu- 
pied in the central region of the up per wire. We note that 
for the experiments of Ref. |24l25j ] , there are at most only 
a few modes occupied in the higher density wire, so the 
assumptions made here may overestimate the quality of 
screening actually obtained in the experiments. 

Let Vo(z,W) — ^-p=5==5 be the interaction with 
short range cutoff W between particles separated by a 
distance z. Here e is the dielectric constant. The en- 
ergy distribution of the wires in Fourier space can be 
estimated as: 



f 



E =- 



dq 



Spu(q)5Pu(q)V (q,W u ) + Spi(q)6pi(q)V (q,Wi) 



2Sp u (q)S Pl (q)V (q,d) 



(24) 



where b~p u /i(q) is the net charge density (electron density 
minus background positive charge) of the upper/lower 
wire in Fourier space and Vo(q,W u /W~i) is the Fourier 
transform of the interaction (already with the short range 
cutoff) within the upper/lower wire. Taking the variation 
with respect to the lower wire charge density 5 pi: 

dE 

-^T=Spi(q)Vo(q,W l )+Sp u (q)V (q,d) . (25) 
Setting it to zero, 



Spi(q) 



$Pu(q)V (q, d) 



(26) 



Vo(<Z, W u ) 

Putting this back into l|24(l and taking out the common 
factors of Sp^, 



E 



S P l(q) 



V (q,W u ) 



v a 2 (q,d) 

V (q,Wi) 



dq . (27) 
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FIG. 3: Effective interaction potential V e a, as a function of 
the separation z between two electrons in the upper wire, tak- 
ing into account the finite width W u of the upper wire and 
screening due to a parallel lower wire of width Wi at a dis- 
tance d. See Fig. Dashed and dotted lines show the bare 
Coulomb interaction Z , and the unscreened interaction soft- 
ened at short distances, l/\/z' 2 + W 2 . Parameters are d = 1, 
W u =0.1 and Wi — 0.15. Region < z < 0.5 is expanded in 
the inset. In the main plot 1/z and 1/V z' 2 + W 2 are indistin- 
guishable in the range 0.5 < z < 4. In the inset the interaction 
potential V c s(z) and \/\Jz 2 + W 2 are indistinguishable in the 
range < z < 0.5. 



Comparing this with the form of self interaction in the 
upper wire it is easy to see the effective form of interac- 
tion in Fourier space is 



V eS (q) = V Q (q,W u ) 



Vg 2 (q,d) 



(28) 



where we have made explicit reference to the short dis- 
tance cut off in real space: Vo(q, W) = dx /^+w 2 ' 
The interaction (|2*5|l is plotted in real space in Fig|3j 



B. Model for confinement and background charge 
distribution 



We now consider the model of a few interacting elec- 
trons in a box of length L. A uniform positive back- 
ground with density equal to the average electron den- 
sity is added to maintain charge neutrality. The resulting 
one-body potential U{x) takes the following form: 



U(x) 



oo, 



-n J Q V c s(x - x')dx', 



if x < or x > L, 
if < x < L. 



(29) 

Here n is the density of positive background charge. 
When the electron number N is fixed, fi = N/L. When 
the electron density changes, as in the case of tunneling 
into a relatively rigid box, n is taken to be the average of 
the initial and final electron density. The precise value 
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of n does not qualitatively change the results discussed 
below. 

Systems with up to four electrons, both with and with- 
out spin, were solved by directly diagonalizing the exact 
many-body Hamiltonian using the Lanczos method^^i 
Earlier applications of the Lanczos method to one- 
dimensional quantum dots were reported in Rcfs. 
liol ] , by Hausler and Kramer, and collaborators. Energy 
spectra for up to N = 4 particles, with Coulomb inter- 
actions unscreened at large distances, were presented in 
Ref. [33, and the splitting between the lowest energy 
singlet and triplet states was found to fall off rapidly 
with increasing wire length in this case. Density distri- 
butions for up to four particles, and density correlation 
functions for three particles, were obtained in Ref. |39| . 
Matrix elements for tunneling were discussed, and tun- 
neling rates were computed, in Ref |4jj > but only for the 
incoherent case, with no momentum conservation. The 
parameters of the calculation, including the interaction 
strength and the geometry of two parallel wires with one 
screening the other, were obtained from estimates in the 
experiments by Steinberg et a/.. 24,25 In the calculation 
shown here we take the upper wire width W u to be 20 
nm, the lower wire width W\ to be 30 nm and the center 
to center distance d between wires to be 31 nm. The 
interaction strength between electrons is obtained by us- 
ing the dielectric constant of GaAs, e = 13.1, and ef- 
fective mass m* — 0.067m e , which gives a Bohr radius 
a B = sh 2 /(m*e 2 ) ps 10 nm. 



C. Ground state density and density-density 
correlation 

Before discussing the momentum structure of the tun- 
neling matrix element, a discussion of ground state prop- 
erties is useful. In the case of electrons with spin, as 
discussed before, the ground state has total spin 0, for a 
system with an even number of electrons, and total spin 
1/2, for a system with an odd number of electrons. As the 
electron density decreases, there is a relative increase in 
the importance of the potential energy compared to the 
kinetic energy, and there is a crossover from a relatively 
uniform liquid state into a Wigner crystal-like state with 
spatially localized electrons. This crossover can be seen 
in Figs. 0] and 01 where we plot the ground state density 
distribution and the density-density correlation function 
for a system of N = 4 electrons. 

In Fig.QJ at a large physical density of twenty-one elec- 
trons per micron, the scaled electron density distribution 
has smooth variation with spatial frequency 2/cp, which 
is simply a Friedel oscillation resulting from boundary 
effects. At the lowest density of one electron per micron, 
the scaled electron density has distinct peaks, whose 
number equals that of the electrons in the system, corre- 
sponding to a Akp Wigner crystal density variation. The 
signature of Wigner crystal correlations is also clearly 
seen in Fig. in the form of 4fcp oscillations in the 
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FIG. 4: Top: Ground state density of a four electron sys- 
tem with spin. Scaled densities are measured in units of 1/L. 
In the label of the individual plot, "density" means the av- 
erage physical density N/L of the system, in units of /jm~ . 
Bottom: Ground state density-density correlation of a four 
electron system with spin. The density-density correlation 
function is defined as ^— ^ J x p(x')p(x' + x)dx' . The pre- 
factor before the integral is to take into account the finite 
length of the system. 



density-density correlation function at low physical den- 
sity. 

The behavior of spinless electrons and spin polarized 
electrons is significantly different from the case with spin 
(shown in Fig. in its spin-singlet ground state. As 
shown in Fig. [5J contrary to the case of a system with 
spin, decreasing the physical density does not signifi- 
cantly alter the ground state. 

Neither the scaled electron density distribution nor the 
density-density correlation function changes more than 
ten percent when the physical density changes over an 
order of magnitude from twenty-one to one electron per 
micron. This means that the spinless system is Wigner 
crystal-like, although essentially non-interacting for all 
densities and gate potentials. This is the result of a com- 
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FIG. 5: Top: Ground state density of a four electron system 
without spin. Bottom: Ground state density-density correla- 
tion of a four electron system without spin. See Fig. 2] for a 
definition of the density-density correlation function. 



bination of the Pauli Exclusion Principle, the cutoff of the 
interaction at large distances due to the second perfectly 
conducting (lower) wire, and the actual parameters of 
the experimental system. When the density is high and 
inter-particle distance z <C as, the interaction energy is 
much smaller than kinetic energy. When the density is 
low and average inter-particle distance z 3> d, the inter- 
action becomes effectively short-ranged due to screening, 
diminishing its effect because in the spinless model the 
wavefunction vanishes when electrons come together. 

As mentioned before, in the experiments of Stein- 
berg et al.M^m as ~ 10 nm and d w 30 nm. Define 
r s = cl/ole, where a is the average inter-particle distance. 
Then, for r s > 3, the long range cutoff set by screen- 
ing from the second wire kicks in and starts to strongly 
suppress interaction effects. Our numerical calculations 
confirm this picture. Since long range interactions are im- 
portant in establishing Wigner crystal-like correlations, 33 
screening from the lower wire acts to inhibit somewhat 



Wigner crystal-like correlations. Longer range interac- 
tions tend to make the peaks in the ground state density 
distribution and in the density-density correlation more 
prominent. Although the effect of the decrease in phys- 
ical density is small in the spinless case, we can see in 
Fig. that distinct peaks in ground state density are 
most prominent in the intermediate physical density of 
eleven electrons per micron. Fig. [5] also shows that at a 
density of eleven electrons per micron, the electrons more 
strongly avoid each other than at both higher and lower 
density. Nevertheless, even at intermediate densities our 
numerical results show that with experimental parame- 
ters of Steinberg et n/.iiiai the interaction effect would 
be quite weak for spin polarized electrons. 

We note that numerical results for the electron den- 
sity and correlation function qualitatively similar to ours 
(inter-wire screening effects not considered) were ob- 
tained in Ref. [4l| . 



D. Momentum dependent tunneling matrix 
element 

Let us now consider tunneling between the interacting 
electron system and a non-interacting electron system, 
shown schematically in Fig. ^ Here we consider the sit- 
uation when the initial and final states of the interacting 
electron system (upper wire) are both ground states. As 
discussed in Sec.[n] the tunneling probability is then pro- 
portional to the matrix element |M(£;)| 2 . 

Typical results of the square of the absolute value of 
tunneling matrix element M are shown in Fig. H3 for elec- 
trons with spin and in Fig. for electrons without spin. 
All results are symmetric under the interchange k — > —k. 
Note that the case with spin shows a much greater change 
with density, although this appears to be limited mostly 
to the amplitude of the peak in the tunneling. 

It is clear from the discussion of previous passages that 
the spinless system ground state under any density can 
be effectively described using a non-interacting model. 
This is clearly seen in the plots. In the spinless model, 
to get the noninteracting matrix element tunneling from 
N — 1 to N electron ground state we put k^ = ttN/L 
in Eq i|15|). Manifestly, this has the consequence that 
M(k = 0) is zero for even N and non-zero for odd N, 
and this is seen in the numerical result. There is also 
a sharp maximum near k^ = ttN/L. As N — > oo, this 
maximum will approach kfj. All these features are seen 
in Figs. The insets of the two plots further demon- 
strate what has already been discussed at the end of the 
previous subsection: Electron interactions cause greatest 
deviation of |M(/c)| 2 from its non-interacting limit at in- 
termediate density. At both very large and small density 
|M(fc)| 2 approaches the non-interacting limit. 

For systems with spin, the interaction has a more sig- 
nificant effect on the matrix element. Yet the salient 
feature here remains that, even at very low density and 
strong interaction, and even after the system changes 
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FIG. 6: Top: Squared matrix elements of tunneling from a 
two to three electron system with spin vs. scaled wavevectors 
kL/n. Here the L = 0.1 /im curve is indistinguishable from 
non-interacting limit. Bottom: Same as Top only tunneling 
from three to four electrons. 

from a Fermi liquid-like state to a Wigner crystal-like 
state, the qualitative features of the scaled matrix ele- 
ment remain similar to those of non-interacting solutions: 
The positions of the maxima are only slightly shifted, 
the positions of zeros a bit more so (apart from k = 0), 
and the values of the maxima are reduced by no more 
than 30% relative to the non-interacting case, for den- 
sities down to 1 [im . Furthermore, since for N = 3 
and N = 4, the systems with spin have the same matrix 
element in the non-interacting limit, the overall features 
of the tunneling matrix elements, after scaling, are very 
similar, even with interactions. 

E. Spin exchange energy 

In the model with spin, the first excited state has to- 
tal spin S = 1 for even numbers of electrons, and while 
for odd numbers of electrons the situation can be more 



FIG. 7: Top:Squared matrix elements for tunneling from a 
two to three electron spinless system vs. scaled wavevectors 
kL/n. Bottom: Same as Top only tunneling from three to 
four electrons. 

complicated.^ 2 The energy gap between the ground state 
and the first excited state decreases as the electron den- 
sity tends to zero. The high density, low interaction part 
of this change in the gap can be understood by simple 
first order perturbation theory using a non-interacting 
wavefunction, as shown in Fig. |SJ In the low density, 
strongly interacting Wigner crystal-like state, electrons 
are spatially confined around their equilibrium positions. 
The system can be approximated as a Heisenberg spin 
chain, where spin exchange is achieved by tunneling be- 
tween adjacent sites iSLSLSi 

N-l 

-ffspin = J ^ Si ■ Si+1 . (30) 
i=l 

This Hamiltonian describes the energy spectrum of low- 
lying spin excitations. Specifically, in the case of N = 2 
we simply have 

J = A , (31) 
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FIG. 8: Comparison of the gap A between the ground state 
and the first excited state energy in a system of two electrons 
in a box, obtained by direct calculation versus by first or- 
der perturbation theory. The unit is Ao, the non- interacting 
energy gap between the ground and the first excited state. 
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FIG. 9: Estimate of the spin exchange energy per electron 
from a two and four electron system using the interaction 
in Eq. I|28ft. Both of the axes are in physical units: the x- 
axis is the physical density, in units of one electron per (im, 
the y-axis is the physical exchange energy J in units of /ieV. 
Typical experimental temperatures are in the range 250 mK-2 
K « 25^eV-200 peV. 

where A is the gap between the ground state and the 
first excited state. For N = 4 numerical solution of the 
four spin Heisenberg chain gives: 

J = 1.51777 A . (32) 

Analytical estimates of J have been obtained in Refs. J3U 
l32H42ll43l | . Hartree-Fock estimates are given in Ref. |30| . 
Related estimates of the spin velocity are given in 
Ref. l44l based on the ladder approximation and in 
Ref. gjj. 

Because the tunneling barrier between adjacent 
sites increases with larger electron-electron interaction 



strength, the effective Heisenberg exchange constant, and 
hence the gap between the ground state and the first 
excited state, will decrease to zero as the physical den- 
sity goes to zero. It is interesting to note that the ef- 
fective Heisenberg exchange energy, extracted from the 
two-electron system and four-electron system gap using 
Eqs. I|31|) and (|32[1 . closely follow each other in the low 
density regime, as seen in Fig. [5] This supports our as- 
sumption that spin excitations in the low-density regime 
are well described by a Heisenberg model, and it sug- 
gests that in this regime, an estimate of the Heisenberg 
exchange from an exact calculation of systems with a few 
electrons can provide a reasonable prediction for the spin 
velocity of an infinite wire with the same electron density. 
It should be noted, however, that the numerical value of 
J may depend sensitively on both short distance and long 
distance cut offsi^ 



F. Mixed spin states 

In the low density regime, where the spin exchange en- 
ergy is small, one may encounter an experimental regime 
where the temperature is larger than J, but still small 
compared to the lowest energy for charge excitations. In 
that case, the tunneling conductance will be a superpo- 
sition of the results obtained from different spin states. 
Also, under experimental conditions, the Zeeman energy 
Ez = g s ptBB may be comparable to or larger than J. 
(Here \xb is the Bohr magneton and g s is the g-factor of 
the electron spin.) If the Zeeman energy is large com- 
pared to both ksT and J, then the electron system will 
be fully spin polarized, and calculations for spinless elec- 
trons will apply. 

As a simplest example, consider the case N — 2, in the 
low-density regime. If the temperature and the singlet- 
triplet splitting A are both small compared to all charge- 
excitation energies in the N — 1 as well as in the N = 2 
states, then the sum in Eq. may be restricted to six 
states: the lowest singlet and triplet states of the N = 2 
system and the two spin states of N = 1 ground state. As 
there are only two independent matrix elements involved, 
the sum may be further simplified to read 

B(k) = w s \M s {k)\ 2 + w t \M t (k)\ 2 , (33) 

where M s (k) and M t (k) are the matrix elements for the 
singlet and triplet states defined in Eq. I|10|) . The weights 
w s and Wt depend on the Zeeman energy Ez and the 
chemical potential fi, as well as on A and T. 

In the limit T — > 0, the ratio w t /w s will approach 
either zero or infinity, depending on whether Ez is smaller 
or larger than A. In the opposite limit, when koF is large 
compared to Ez and to A, we have Wt/w s — > 3/2. In Fig. 
I1UI we show the momentum dependence of B(k) for these 
three limiting values of wt/w s , at L = 0.4 /im. 

To illustrate further, we may consider the case where 
A is small compared to Ez and to fc^T, but the ratio 
Ez/iksT) is finite. If [i is adjusted so that the upper 
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wire has equal probability of having N 
electrons, then we find 



1 or N 



{2x 2 + 5x + 2) 
6x 



(34) 



where x = exp[i?z/(fcs7 1 )]. As expected, we see that 
Wt/w s = 3/2, for x = 1, and w t /w s — > oo, when x — > oo. 



singlet 
triplet 
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FIG. 10: Momentum-dependence of tunneling conductance 
between one-electron and two-electron states for a wire of 
length L=0.4/xm. Dashed curve shows line-shape |Af s (fc)| 2 
for tunneling into the singlet ground state of the two particle 
system, applicable at T = when the Zeeman energy Ez is 
smaller than the exchange gap A. Short-dashed curve shows 
\Mt(k)\ 2 , for tunneling into the triplet ground state, applica- 
ble at T — when Ez > A. Solid curve is a weighted average, 
applicable if T is large compared to both Ez and A but small 
compared to the energy of the lowest charge excitation. 



VI. LIMIT OF LARGE N 

Properties of an infinite one-dimensional electron sys- 
tem have been explored by a variety of techniques, includ- 
ing solutions of exactly solvable models, renormalization 
group methods, bosonization, and conformal field the- 
ory. These techniques can be also adapted to the study 
of finite systems, when N is large. 

In this section, we will use mostly bosonization tech- 
niques to explore both spin-coherent and spin-incoherent 
regimes. Bosonization techniques give most readily the 
electron Green's functions 



g a {x,x',t) = (Vv(M)Vi(z',o)> 



(35) 



as a function of position and time t. For an infinite 
system, position variables enter only in the combination 
(x — x') 1 but for a finite system, the distance from the 
ends may also be important. To obtain the spectral den- 
sity at a finite energy %lo (setting \i = 0), one must take 
the Fourier transform with respect to t; however one may 
generally obtain an estimate of this by evaluating G at 



an imaginary time r equal to uj^ 1 . To obtain the ground- 
state tunneling amplitude for a finite system, when fceT 
is small compared to all excitation energies, we need to 
study the Green's function for r — > oo. In particular, 
if the energies are adjusted so that E N = fi™" 1 , then 
evaluating Ij35(l at a (N — l)-electron ground state yields 



lim g(x,x',T)=^ s (x)*? s *(x') 



(36) 



where fy^(x) 



is the quasi-wavefunction defined in 
Eq. (|12(l . In the free-spin regime, at finite temperatures, 
we shall make use of a generalization of this result. 



A. Ground-state tunneling regime 

Numerous works have considered tunnel- 
ing in the regime where the spin energy is 
important.?! 12 ' 14 ! 1 '?! 21 ' 22 ! 2 ?! 41 ? At sufficiently low en- 
ergies, an infinite interacting one-dimensional electron 
system is believed to be described as a Luttinger liquid^ 
in which spin and charge propagation are characterized 
by two distinct velocities, v s and v c . At finite bias volt- 
ages, the amplitude for momentum resolved tunneling 
into a Luttinger liquid is expected to show structure 
reflecting both of these velocities. This follows from the 
form of the one-particle Green's function (z = x — x') 



G{z,t) 



y/(v s T - iz)(v c r - iz)[(v c t) 2 + z 2 } 



21 n 



dropping a short-distance cutoff-dependent prefactor. In 
the limit where the bias voltage V — > and T — > 0, 
however, the momentum-resolved tunneling conductance 
reduces to 8 functions at k = ±kp similar to that of 
a non-interacting electron system. The primary effect 
of the interactions is then to reduce the amplitude of 
the (5-functions. For an infinite system, the amplitude 
is predicted to vanish as max(/csT, V) a , when T and V 
approach 0, where the bulk tunneling exponent a > de- 
pends on the strength of the electron-electron interaction. 
Tunneling into a point near the end of a semi-infinite sys- 
tem is characterized by a different exponent cv clK j, which 
is generally larger than a. The tunneling exponents are 
related to the Luttinger liquid interaction parameter g 
by a = (g + g~ x — 2)/4 and a cn d = (g^ 1 — l)/2 for a sin- 
gle mode wire^i For fermions with repulsive interactions, 
as considered here, one has g < 1. For translationally- 
invariant systems, g = vf/v c , where vf is the Fermi 
velocity with turned off interactions. For not too strong 
interactions, one can make an RPA estimate for g of a 
cylindrical wire of radius R with electron density n at a 
distance D from a two-dimensional screening gatei^2> 




8\n(2D/R) 

■K 2 na,B 



(38) 
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We remark that if the lower high-density wire had charge 
modes propagating with a very different velocity from the 
upper wire, the interwire interactions would not signifi- 
cantly couple or modify the charge modes of individual 
wires— Using parameters of Ref. 24.25], 2D ss 1 /im, 
R f» 20 nm, as ~ 10 nm, and n « 30 jmr 1 near the lo- 
calization transition, we get g ~ 0.3 which is comparable 
to the measured g ss 0.5. 

Kane, Balents, and Fisher— considered tunneling into 
a finite-length metallic carbon nanotube, which has four 
excitation modes, rather than the single charge and spin 
modes considered here. Nevertheless, their formulas may 
be readily adapted to the present case. The general re- 
sult, for a wire with uniform electron density and hard- 
wall confinement at the ends, may be written in the form 



*cff 



(x) 



Vln~ c 



nx\ 

TJ 



ij(a on d-a) 



sm(k F x) , (39) 



valid for positions not too close to the wall where the 
factors involving a cease to be a good approximation. 

As a increases with increasing interaction strength, 
the overall amplitude of the tunneling matrix element 
M(k), obtained by Fourier transforming (|39|l . will de- 
crease rapidly with increasing interaction strength, for 
N large but fixed. However, the momentum-dependence 
remains qualitatively similar to the form (|15fl for non- 
interacting electrons, and the peak at kp broadens only 
gradually with increasing interaction strength, as illus- 
trated in Fig. ^2 
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FIG. 11: Momentum-dependence of ground-state tunnel- 
ing for N=10 electrons and various interaction strengths, ob- 
tained by calculating |M(fc)| 2 from the effective wavefunction 
Each curve has been divided by a normalization con- 
stant N g , chosen so that the plotted curves have equal areas. 
The norm N g decreases rapidly with increasing interaction 
strength (decreasing g), but the width of the peak broadens 
only slightly. Note that for g — 1 the results are equivalent 
to the N = 5 case in Fig. [5] 

Tserkovnyak et alM^ have explored the implications 
of Luttinger liquid theory for momentum-resolved tun- 



neling into a long finite wire with soft confinement at 
the ends. Although the emphasis of that work was on 
tunneling at finite bias, where the discreteness of the 
energy levels plays a relatively minor role, the analysis 
can also be applied to estimate the form of ty^ s (x) and 
the matrix element M(k) for the ground-state tunnel- 
ing into a wire with large N, for not too strong interac- 
tions. Essentially, the factor sm(kpx) in (|39|l should be 
replaced by the WKB wavefunction for a non-interacting 
electron in a self-consistent potential that gives the cor- 
rect density profile along the wire. The zeroes of ^f^ s (x) 
will be spaced further apart near the ends of the wire 
than near the center, but near the center of the wire, the 
wavefunctions will look quite similar for the cases of soft 
and hard confinement. For increasing values of the the 
electron-electron interaction, the envelope factor in l|39f) 
is increasingly weighted towards the center of the wire, 
so that the Fourier transform M(k) should become corre- 
spondingly insensitive to the difference between soft and 
hard confinement. 



B. Finite temperatures 

At finite temperatures T, for the infinite system, one 
expects that the 8- functions at k = kp will be broad- 
ened by an amount of order 8k rj ksT/ (hv s ). If the 
temperature is larger than the energy of the lowest spin 
wave mode, irhv s /L, the broadening will be larger than 
2ir/L for a finite system of length L. The width 8k be- 
comes comparable to kp when fc^T becomes comparable 
to J ~ hv s kp 7 the energy of the shortest wavelength spin 
excitations. 

If the spin velocity is very much smaller than the 
charge velocity, one can enter a regime where fc^T is 
larger than the energy of the shortest-wavelength spin 
excitations, but still low enough that relatively few 
short-wavelength charge modes are excited. This regime 
has been explored in recent papers by Cheianov and 
Zvonarevf^Si and by Fiete and Balents— The form fac- 
tor for momentum-resolved tunneling is indeed broad- 
ened by an amount of order kp in this regime, as the 
Green's function in position space is found to fall off ex- 
ponentially with distance. A more precise description of 
this behavior will be given in the next subsection. 

The lowest energy cost for a charge excitation is e c = 
nhv c /L. The highest energy for a spin state with mul- 
tiple excitations is ~ Nhv s kp. Earlier, we defined an 
extreme free-spin regime as having a temperature high 
enough so that all spin states have equal weight, while 
there are no charge excitations. This makes sense for a 
small system, but is very restrictive for a large system, as 
it requires that v s /v c <C l/N 2 . However, the spins will 
actually be effectively free as long as fcgT is greater than 
the energy of a single short wavelength spin-excitation, 
while smaller than the lowest charge excitation. This 
only requires that v s /v c <C l/N . The results for the mo- 
mentum dependence of the tunneling will be the same in 
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this case as for the extreme free-spin regime. 

C. Free-spin regime 

We first focus on the regime with energy-scale hier- 
archy J <C ksT <C hvcttF, so that all spin configura- 
tions have effectively equal statistical weight while short- 
wavelength charge excitations are absent in equilibrium. 
(We do not yet impose the more restrictive condition 
kgT <C e c .) We are particularly interested in the Green's 
function Q(x,x',t) at large imaginary times r but spa- 
tial separations of the order of the inter-electron spacing, 
as these are the distances that dominate the tunneling 
response, as discussed in the following. For times smaller 
than the time for a charge excitation to cross the system, 
G should depend only on the separation x — x'. For longer 
times, there may also be a dependence on the distance 
from the boundary. We will first calculate the Green's 
function in the infinite wire limit. 

A quantity that will play a crucial role in our analysis 
of the tunneling at r 3> a/v c is the root-mean-square dis- 
placement of an electron during a time r due to quantum 
fluctuations [derived in Eq. I]54[l]: 



u(t) = - ^/2g\n(v c T/a) , (40) 

where a is the mean spacing of the electrons. When 
u(t) <C a (which is possible only for very strong inter- 
actions, g -C 1), electrons can be thought of as forming 
a "true Wigner crystal" on the relevant energy scales. 
When u(t) > a (which for strong interactions corre- 
sponds to exponentially small energies), fluctuations de- 
stroy the long-range translational order, but significant 
Wigner crystal correlations may be present at intermedi- 
ate length scales. We refer to this situation as the "fluc- 
tuating Wigner crystal" . For a finite size system, the 
fluctuations u(r) are cut-off by the length of the system 
at times r = L/v c . The momentum structure of the tun- 
neling will exhibit qualitative differences depending on 
the relative size of a and u, and the form of the Green's 
function itself will also depend on the relative sizes of a, 
u, and x — x' . 

The condition J <C ksT <C hv c kp requires that 
the short-distance repulsion between electrons be very 
stron g 31 i , so that neighboring electrons exchange po- 
sitions only on a time scale h/J which is longer than 
t ~ h/(ksT). The Green's function Q can then be ex- 
pressed as a series of the form: 

g(x,x',r) = J2Sjm(x,x',T)(-l) m Pm , (41) 

jm 

where Gjm is the Green's function for spinless bosons 
which do not change places, when a particle is inserted at 
time at space point x' as the j th particle on the lattice, 
and removed at time r at point x in lattice position j+m. 
The factor (— l) m encodes the Fermi statistics associated 



with the permutation of a given electron through m inter- 
mediate electrons, and if we were to set p m = 1, we would 
recover the Green's function for spinless fcrmions. Here, 
however we set p m — 2~' m ', which is the probability of 
finding \m\ + 1 consecutive spins pointing in the same 
direction. In the limit we are considering, where all spin 
states are effectively degenerate, and the spin dynamics 
is effectively frozen out (exchange events do not occur), 
all electrons between lattice positions j and j + m must 
have parallel spins to contribute to the Green's function 
(i.e., for the world line trajectories to wrap around the 
imaginary time torus, as discussed in Ref. 51]). 

Let us first consider the situation where \x— x'\ is much 
larger than a and u. The creation operator for a spinless 
fermion can be bosonized in terms of 9 and <fi fields asso- 
ciated with the fluctuations in the electron density and 
momentum density. 

Specifically, 9 is related to the particle density fluctua- 
tions through the familar relationship 6n(x) = -^d x 9(x), 
and the fields satisfy [cj)(x),9(x')) = —mlKx — x) where 
H is the Heaviside function. In a first-quantized path 
integral formulation of the Green's function, the domi- 
nant contributions come from those world-line trajecto- 
ries that fluctuate only a small amount from configura- 
tions corresponding to the mean density (due to strong 
Coulomb repulsion in the Wigner crystal limit), as dis- 
cussed in Ref. |51| . 

For a, u <C \x — x'\ <C L and x, x' away from the wire 
boundary, density fluctations can be treated as a con- 
tinous variable and the Green's function (|41|l may be 
expressed as^i 

Q(x x' t) ~ (2~\ n ( x > x '' tS >\(— i^ N ( x ' x '> T ) e l l<P(. x > T )-<P( x 'M l j 

(42) 

Here N(x, x' ,t) is the particle number operator between 
points x and x'. The first factor, in (|42[) corresponds 
to p m , while the second and third together represent the 
product of creation and annihilation operators for a spin- 
less fermion. The average is taken at T — > in the charge 
sector. Using the bosonized expression of the particle 
number in a distance x — x' > 0, 

N(x, x', t) = n(x - x') + - [0(x, t) - 6(x' , 0)] , (43) 

the electron number can be written in terms of the aver- 
age density, n, and a fluctuating piece expressed in terms 
of the 9 fields. Using Eq. Ij43(l . and approximating op- 
erator (— 1) N as Re[e i7rAr ], the Green's function can be 
expressed as Q = Q+ + Q- , where^i 

g + (x,x',T) ^e-'H^^'lve^^^Oie-irPM- 9 ! 1 '")] 

x e i[e(x,r)-e(x',0)] e i[^(x,T)-<j>(x',0)]\ ; (44) 

and t/„ (as we will see from the derivation below) is given 
by Q- = [{?+]*. The first two factors in Eq. g2J come 
from the exponentiation of the average density and we 
have defined fc^ = nh. The exponential decay of the 
first factor results from spin averaging and the oscillatory 



14 



second term comes from Fermi statistics, as discussed 
above. 

We now compute the part of the Green's function com- 
ing from fluctuations in the charge sector using the effec- 
tive low energy "elastic" action 



S r = I dxdr— 
2?r 



^{d x 9f + 2g(d x 



-d T c, 

IT 



obtained from the Hamiltonian 



H = \ dx— 

2tt 



±(d x 8) 2 + 2g(d x 0) 2 
2 9 



(45) 



(46) 



A justification of the Harmonic form 146JI for a stro ngly 
interacting system has been given in Refs. [23 and |52j |. 
Since the action is quadratic, the average can be moved 
to the exponent 



G+(x,x',t) 



-k F (x-x')i^ ik F (x-x') -^(l+i^) 2 (0 2 ) 



x e 2 



l<* 2 )p-(i+^X*e> 



(47) 



where $ = <p(x,r) - (f)(x',0) and 9 = 6(x,t) - 0(x',O). 
If we assume as before that the electrons are confined to 
a wire of length L, the fields must satisfy the boundary 
conditions d x <f>{x = 0, L) = and 6(L) - 0(0) = n(N - 
nL). Here iV is the total number of electrons. Expanding 
the fields in a Fourier series gives^ 



OO 

E 

m— 1 



2g . fm/n.v \ 
i\l — sin 

m V L 



) Oh, 



bl 



1 



-e^(x) (48) 



Green's function Q = Q + + Q_ was computed earlier i^SiSi 



G(z,t) 



C'e-' k ^^ 
(z 2 + v?t 2 ) a o 



i(k F z~tp+) 



V C T — IZ 



-i(k F z-ip g ) 
V C T + iz 



(51) 

where C is an undetermined constant. The phases <p~ 
are given by 



<Pg ( Z , T ) 



In 2 

7T 



.9 l n 



Z 2 + V 2 T 2 



1 , / v c t — iz 
± - In — 

2 \V C T + IZ 



(52) 

The power law decay of (|51(l is characterized by the 
anomalous exponent 



1 

89 



1 - 



In 2 

7T 



(53) 



We note that the Green's function H51fl is valid only for 
z a,u and that it decays very rapidly with distance 
due to the incoherent spin degrees of freedom. The ex- 
ponential piece in (|51(l shows that the spin incoherent 
regime (J <C fc^T) destroys the strictly power law decay 
of the correlation function present in the Luttinger liquid 
theory where the spin degrees of freedom are coherent. 
In order to calculate the momentum resolved tunneling, 
we shall need an alternate formula for C/(z,r), which re- 
duces to (|51|l when a, u <C z <C v c t, but is also valid 
when z ~ a, u. 

When the distance z is comparable to the separation 
between electrons (where the density fluctuations are 
large compared to the average density), but z <C v c r 
we have in the limit L — > oo^i 



(x,t) 



E 



V 2gm 



1 /rmrx\ ,, 

cos ( — ^- J (b n 



L 



bl 



+$ c , (49) 



where the zero mode term 6^ \x) = j;ir(N — nL), the 

b m satisfy \b m ,b^ m ,] = 6 mm > and the operators N and <E> C 
satisfy [N, $ c ] = 1. These field expansions diagonalizc 
the Hamiltonian (|46() of the charge sector as 



00 



H = ^-N 2 + me c blb m 

rn — 1 



(50) 



The first term is the ground-state energy in the charge 
sector, with the charging energy Ec = e c /(2g). The sec- 
ond term represents the energy associated with collective 
density oscillations of the low density electron gas. 

Continuing with the evaluation of the Green's function 
for z — x — x' 3> a, u, it is useful to first consider the 
limit of L — *• 00. Then the correlators in the exponent of 
(|T7jl depend only on z and can be evaluated as (6 2 ) = 
gln[(z 2 + v 2 T 2 )/a 2 }, ($ 2 ) = (Ag)- 1 \n[{z 2 + v 2 t 2 ) / a 2 } and 
($9) = 2~ 1 ln[(v c r - iz)/(v c T + iz)]. (Here a = n' 1 
is the average spacing between electrons.) The resulting 



G(z,t) 



E 

j=-oo 



2 _lil {-l) j (5{N{z, t) - jyW--,r)-^o,o)\) 



OO 



^2 _ IJ'I(— iy J — e - iX (3-nz)^ e i[Xe(z,T)/7T+i(z,T)]^ 
j=-oo 

= ^-b-|(__l)y^ e -^o--fi Z ) e -^(^<e 2 >+2A(e*>+<* 2 

(54) 



— y 

j=-co 



o 2—(-^y 9 f(z,u(T)) , 

2ttu(t) \v c t ) 



where 



f(z,u) 



00 

E 

J=-oo 



2-1^1 (-l)J'e" 



(55) 



In the second to the last line of (|54l) , we have made use of 
the condition z <C v c t to neglect the contribution of the 
correlator (9$), and in the last line we have evaluated 
the remaining correlators and used the definition (|4UII for 
u. 
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With Eqs. (|51(l and (|54(l in hand, we are now ready 
to study the momentum-resolved tunneling spectrum 
of a long quantum wire in the limit J <C fc^T <C 
hv c kp, by computing the Fourier transform Q{k,r) = 
J dze~ lkz Q(z,T). We are interested in large values of t, 
so the Fourier transform will be dominated by contri- 
butions from z <C v c t. The dominant, short-distance, 
z ~ a, correlations that determine the momentum com- 
position of the Green's function should thus be correctly 
described by Eq. (|54|l . The Fourier transform of Eq. Ij55(l 
can be computed directly, but it is convenient to express 
/ as the convolution of a Gaussian with the function 



(56) 



j = -oo 



The Fourier transform of / is thus the product of a Gaus- 
sian and the Fourier transform of l(z), which leads to 



G(k,T)~a[ — I'" A(M(r)) > (57) 

1 Vr.T ' 



A(k,u) 



5 -I- 4cos(fca) 



(58) 




FIG. 12: Spectral function A(k, u), which determines the mo- 
mentum dependence of tunneling in the free-spin regime. The 
quantity u is the root-mean-square electron displacement, due 
to quantum fluctuations, from the sites of a classical Wigner 
crystal, and a is the lattice spacing. When u > a, the mo- 
mentum distribution is single lobed and peaked about zero 
momentum. In the opposite limit, when u a, the mo- 
mentum distribution exhibits a doubled lobed structure with 
peaks near k — ±fci? = ±7r/a. The peaks have width slightly 
less than k F , in agreement with Eq. 16011 obtained from the 
Green's function 15111 . 



The amplitude for momentum-resolved tunneling into 
an infinite wire, at an energy %u>, is determined by the 
inverse Laplace transform of the function Q(k,r). In 
practice, this is determined by the form of Q{k,r) for 
r m lu^ 1 . For tunneling into the charge ground state, we 
shall actually want to take the limit r — * oo, with the 
system length L held finite. We shall see below, that the 
behavior of Q in this limit can be obtained, to a good 
approximation, by using the expressions for an infinite 
system, with r replaced by L/v c . Thus the momentum 
dependence will be given by the function A{k, u) defined 
in H58|) . with u evaluated as 



u(L/v c ) = -y/2g\n(L/a) 



(59) 



Equations (|57|l and (|58|l are, therefore, central results 
of this section and they have several features worth em- 
phasizing. The first is the momentum structure: There 
is an exponential envelope centered about zero momen- 

k 2 a 2 

turn, e ~ , whose width is given by the parameter 
u measuring the fluctuations of an electron's position. 
Larger fluctuations imply a more sharply peaked enve- 
lope in momentum space. This envelope multiplies an- 
other momentum dependent function, which is sensitive 
to the mean spacing of the electrons and has maxima at 
k = ±7r/a = ±fcp. 

Clearly, the overall shape of the momentum distribu- 
tion will depend on the relative size of u and a. The 
maxima at finite k will disappear and merge into shoul- 
ders for u/a larger than a critical value of order 0.75. By 



contrast, in the limit u/a —> 0, the height of the max- 
ima is nine times the value at k = 0. The momentum- 
dependence of the spectral function A(k 7 u) is shown in 
Fig. E| for several choices of u. 

One can also obtain the double-lobed structure of 
G(k,r) in the limit where u <C a by Fourier transform- 
ing Eq. Q5ip. although it is, strictly speaking, only a 
valid approximation when z 3> a, while it is the re- 
gion z ~ a which dominates the spectral properties. 
Neglecting the slowly- varying phases y£ , Eq. (|51l) gives 
G(k,r) w A+ + A-, with 



A±{k) oc 



1 



ln2\2 



(k±k F ) 2 + (k F ^ 



(60) 



The density of states for tunneling into a single point 
x, far from the ends of the wire, at a finite energy ui, 
is approximately given by A{uj) ~ tQ(z = 0, r), with 
t = lu^ 1 . Using (|54|l . we see that when u is of order 
a or larger, the frequency-dependence of the tunneling 
density of states is given by A(u>) oc oj a /u, with^i 



1 

a=—-l. 
4.9 



(61) 



This implies that for Ag > 1 the density of states di- 
verges as u> decreases, unlike the behavior of a Luttinger 
Liquid, while for 4g < 1 the tunneling density of states 
is suppressed, qualitatively similar to the behavior of a 
Luttinger Liquid^ 

The density of states must be modified if one tunnels 
into a point near to the end of a semi-infinite wire, so that 
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the distance from the end is small compared to v c Ju>. For 
r > x/v c , the expectation values (6 2 ), ($ 2 ), and ($6), 
which enter l(54|) are affected by the proximity of the 
boundary. In particular, fluctuations in are reduced, 
but ($ 2 ) is doubled near the end of the wire, leading to 
a reduction of the tunneling amplitude. One finds, for 
fixed value of x and to — > 0, that the tunneling density of 
states scales as A(x, uo) cx u a =»<ij; a ™4- a ) where 



<5end = 7. 1. 

25 



(62) 



We now discuss in greater detail the effects of finite 
system length L. To obtain the resonant tunneling am- 
plitude in the free-spin regime, we need to calculate the 
function 



g(x,x') = [J2(P\Mx)\<x){<x\tiWM} 



(63) 



Here, a and (3 denote the spin states of the system with 
N and N — 1 electrons, respectively, and the average is 
taken over the spin states j3. The tunneling conductance 
(Q) is given by [see Eq. (JSJ] 



\ 2 Pv / dxdx'e- lk±< > x - x ")g(x,x') 



(64) 



The function Q(x 1 x') is determined by the function 
Q(x 1 x' ,t), in the limit r — + 00, which we can evaluate 
using the same steps as in Eq. (15411 above. Because of the 
exponential dependence of the first factor in Eq. I|47|) , we 
are only concerned with situations with \x — x'\ <C L. 
We must now take into account the discreteness of the 
normal modes, as well as the proximity of the boundary, 
with the result that the expectation values (0 2 ) and ($ 2 ), 
entering 154(1 . depend on xj L, as well as, logarithmically, 
on L/a. The general result is 



G{x,x') 



L- 1 {a/L) & 
\J9 ln(i/a) 



sin 



(f) 



f(x-x'), 



(65) 

where a and a on d are given by (|61J) and l(62|l . and f(x—x') 
is given by l|55|) with u = u(L/v c ) . (The derivation of 
these results is similar to that used in Ref. to discuss 
tunneling into a finite carbon nanotube.) 

The amplitude for momentum-resolved tunneling is de- 
termined by the Fourier transform of l|65l) with respect 
to the difference variable (x ~ x'). Thus we recover the 
/c-dependence A(k,u), given by with u = u(L/v c ). 

We see that in the limit of large N there are two qual- 
itatively different features to be expected in the tunnel- 
ing between parallel quantum wires in the free-spin (spin 
incoherent) regime compared to the ground state (spin 
coherent) tunneling regime: (1) The momentum struc- 
ture is no longer sharply peaked at k = ±kp at zero 
bias with the peaks splitting at the spin/charge velocity 
slopes at a finite bias, but depends on the relative size 
of u and the mean electron spacing a. When u»a, the 
momentum distribution has a single broad peak centered 



at k = 0, while when u <C a, the momentum distribu- 
tion has a double lobed structure peaked at k = ±kp at 
zero bias with the peak position shifting at the charge- 
velocity slope at a finite bias. Note, also, that kF is twice 
as large as /cf for unpolarized electrons at the same den- 
sity. (2) Depending on the interaction parameter g of the 
charge sector, it is possible to have a diverging tunneling 
density of states as the energy is lowered. This behavior 
contrasts with the ubiquitous power law suppression of 
the spin-coherent Luttinger Liquid regime. 



D. Effects of non-zero Zeeman field 

Our results for the free-spin regime can be generalized 
to the situation where the Zeeman energy Ez = g s fJ-BB 
is comparable to ksT 1 but large compared to the ex- 
change energy J. In this case the spins will be partially 
polarized, but uncorrelated, and we can use essentially 
the same analysis as before. The primary modification 
is that the probability p m for |m| + 1 consecutive spins 
that are all spin up or all spin down, respectively, which 
appears in Eq. (|41() . is now given by 



P« 



w 



M+i 



\m\ + l 



(66) 



where w+ = 1— w_ = (l+e~ Ez / fc BT-pi i s the probability 
for finding a spin aligned parallel to the field. 

To implement this change, the factors 2"' in (|54() and 
(|55|l must be replaced by pj . The function A(k, u), which 
describes the momentum dependence of the tunneling 
rate, is now given by 



A(k,u) 



F 4 



72 



«4) 



2w± cos(ka) 



(67) 



(68) 



It may be readily seen that when the degree of polar- 
ization is increased, the peaks of A(k,u) at k = ±kp 
become higher and narrower, while the weight decreases 
elsewhere. (The higher momentum peaks are exponen- 
tially supressed by the prefactor in Eq. I|67| l.1 If ka is 
fixed at any value other than an odd multiple of 7r, then 
the value of A(k,u) will vanish in the limit w + — > 1. 
By contrast, if ka is equal to an odd multiple of it, the 
value of A{k, u) will diverge in this limit, proportional 
to wZ 1 ■ The function A(k,u) is plotted in Fig. 1131 for 
several values of x — Ez/ksT with u — 0.6a. Compare 
with Fig. H2 

To estimate the amplitude for resonant tunneling into 
a single point near the center of the wire, we can use l|65|l . 
with x = x' = L/2 . The function / is to be evaluated 
using l|55(l . with the factor 2^^ replaced by pj, and u 
given by (|59|l . For any value of the polarization less 
than unity, the value of f(z = 0,u) is finite in the limit 
u — > 00 , so we find 



Q(x, x) — L 



-(1+5) 
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-1/2 



(L/a) 



(69) 
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FIG. 13: Polarization dependence of the spectral function 
given by Eqs. and JfjSJ. Here x = E z /k B T with E z 

the Zeeman energy and T the temperature. We have used 
u — 0.6a for each curve. (See Fig. 1121 for the dependence on 
u at Ez = 0.) As the polarization increases, the peaks at 
fci? = n I a become narrower and higher as described in the 
text and the weight between the peaks decreases. 

for L — > oo with a given by l|61|l . On the other hand, 
for a fully polarized system, with w + = 1 , one finds that 
f(z = 0,u) vanishes like uL~~ 9 for L — ► oo and u given 
by One then finds 

G(x,x) ~ 1/^+^) , (70) 

with a p = g + (4g) _1 — 1 . This exponent is the same 
as the standard tunneling exponent one would obtain for 
spinless fermions having the same interaction and den- 
sity, if one takes into account the difference in the defini- 
tions of g for the two cases. (Here, we have defined g as 
vf/v c , where Dp, the Fermi velocity for non-interacting 
spinless fermions, is half of the Fermi velocity for spinless 
fermions at the same density.) 

VII. COMPARISON WITH EXPERIMENTS 

Results for the momentum-dependence of the tunnel- 
ing conductance, obtained in the present paper, show 
various similarities and differences with the experimen- 
tal results of Steinberg et alM^ in the localized regime 
(that is, the regime in which the tunneling exhibits broad 
structure in momentum space and hence local structure 
in position space). In general, experimental results at 
the first Coulomb blockade peak, associated with the first 
electron to tunnel into the localized region, show a single 
peak, centered at k = 0, as expected from our analysis. 
However, the spectrum for the second and subsequent 
electrons typically show two peaks, at positions rs ±fc m 
that increase with increasing N. The intensity between 
the two peaks is not small, being perhaps half the inten- 
sity at the maxima, and one does not observe the zeroes 



of the intensity that we find for ground-state tunneling 
in a symmetric well. 

Qualitatively, the experiments are most consistent 
with our results for the free-spin regime, i.e., the results 
one would expect when fc^T is large compared to J but 
small compared to the energy for charge excitations. Ac- 
cording to the results of Subsection IV El above. J should 
indeed be smaller than the experimental temperature of 
0.25 K, for densities less than about 6 electrons per /jm. 
The confinement length L in the experiments is not pre- 
cisely known, and it may well change as additional elec- 
trons are added, but if one estimates L in the range of 
0.5 to 1 /ion, the free-spin regime would hold for values 
of N less than perhaps 3 to 6. On the other hand, there 
is evidence that for the largest values of N in the local- 
ized regime, the density reaches 15 electrons per /im, for 
which we estimate a value of J much larger than fc^T. 

An added factor in the experiments is that the Zeeman 
energy is not negligible compared to fc^T. For a typical 
applied field of 2 T, the Zeeman energy is about 0.6 K, 
more than twice the quoted temperature. If the Zeeman 
energy is comparable to, but smaller than J, the system 
will be partially polarized even at T — 0. For J ~ Ez > 
ksT > J/N, typically a number of different spin states 
will be occupied. Theoretical expectations for the case 
with a sizable Zeeman energy were discussed briefly in 
Subsection lVI Dl above. for the free spin regime, and they 
may be qualitatively consistent with the experiments. 

We recall that in the free-spin regime (provided that 
u < 0.75a), or in the fully-polarized regime, the position 
of the intensity maximum should occur at k m « 2/cp, 
where kp is the Fermi wavevector for unpolarized elec- 
trons at the same density. There is no experimental in- 
dication of this when the upper wire is in the delocalized 
(that is, the Luttinger Liquid regime) regime. Thus if 
the electrons in the localized regime (see first paragraph 
of this section above) form either a free-spin or a spin- 
polarized lattice, there should be an increase in the value 
of k m as one enters the localized regime. It is not clear 
whether such an effect is observed in the experiments. 
The tendency for an increase in k m may be counteracted 
if there is a discontinuous decrease in electron density 
upon entering the localized regime. 

The situation where fc^T and Ez are both large com- 
pared to J was discussed for the simplest case, N = 2, 
in Subsection IV Fl and illustrated in Fig. ^| Under 
the experimental conditions where Ez is large compared 
to ksT, the momentum dependence of the conductiv- 
ity should be close to that of the pure triplet state, and 
should have near zero intensity at k — 0, if the confining 
well is symmetric. This contrasts with the experimental 
data for N = 2, where the intensity at k = is a substan- 
tial fraction of the maximum intensity. The discrepancy 
might be explained if the confining well is sufficiently 
asymmetric, or if the electron temperature in the upper 
wire is significantly higher than the lattice temperature. 

The calculations in Section V above used a model with 
sharp confinement at the ends of the wire and a flat po- 
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tential within. We presented exact numerical diagonal- 
izations for two and four electrons. A model with soft 
confinement, and higher density near the center of the 
wire, would lead to more intensity at \k\ < k m , and a 
faster fall-off of intensity for \k\ > k m , than one obtains 
in the case of sharp confinement 

Finally, we repeat that the reason for the apparent for- 
mation of a barrier between the low-density central sec- 
tion of the upper wire and the high density outer regions 
is not well understood. 



VIII. CONCLUSIONS 

In this paper we have studied momentum resolved tun- 
neling into a short quantum wire. We have focused on 
the situation where the number of electrons, N, is not 
too large and the density of electrons is low enough to 
approach the Wigner crystal limit. We began by review- 
ing some general theorems which dictate basic properties 
of the momentum dependence. We note that the matrix 
element for tunneling between ground states for N and 
N — 1 electrons can be expressed in terms of the Fourier 
transform of a real function 'Jeff (a;), which is the ma- 
trix element of ip a (x) , the electron annihilation operator 
at point x, between the two states. For non-interacting 
electrons, i& e g(x), is the wavefunction of the highest filled 
level in the iV-electron system; more generally, fy e g(x) 
is a quasi-wavefunction, whose magnitude is reduced by 
many-body correlations, but whose spatial dependence is 
qualitatively similar to the non-interacting case. 

Using exact diagonalization, we computed the ground- 



state density, the density-density correlation function, 
and the effective exchange constant, J, for N < 4. We 
include screening effects from the adjacent, higher den- 
sity wire. In the limit of large N wc applied bosoniza- 
tion techniques to study the momentum resolved tunnel- 
ing, both at T = 0, and in the free-spin limit, of spin 
energies much less than temperature and the Fermi en- 
ergy, J < k B T < Ep. Whereas for k B T < J,E F , 
the momentum dependence has narrow peaks centered 
at k = ±kp, in the free-spin regime, the spectral func- 
tion can have broad peaks centered at k — ±2fcp. Com- 
parisons of our theoretical predictions with experimental 
results obtained by Steinberg et ah 24 i 25 show qualitative 
agreement on many features, but some aspects are still 
not understood, and further work is needed. 

While we were finishing this manuscript, we became 
aware of independent unpublished work by E. J. Mueller, 
based on the experiments that motivated our study. 54 
Mueller's work employs, primarily, a Hartree-Fock ap- 
proximation, and is largely orthogonal to the work re- 
ported here. 
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